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ABSTRACT 

Templates for polarised emission from Galactic foregrounds at frequencies relevant 
to Cosmic Microwave Background (CMB) polarisation experiments are obtained by 
modellin g the Galactic Magn etic Field (GMF) on large scales. This work extends the 



results of |Q'Dea et al. ( |2Q12[ ) (hereafter FGPolI) by including polarised synchrotron 
radiation as a source of foreground emission. The polarisation direction and fraction in 
this calculation are based solely on the underlying choice of GMF model and therefore 
provide an independent prediction for the polarisation signal on large scales. Tem- 
plates of polarised foregrounds may be of use when forecasting effective experimental 
sensitivity. In turn, as measurements of the CMB polarisation over large fractions of 
the sky become routine, this model will allow for the data to constrain parameters in 
the, as yet, not well understood form of the GMF. Template foreground maps at a 
range of frequencies can be downloaded from the on-line repository [f] 

Key words: cosmic microwave background, polarisation experiments, foregrounds, 
£?-modes, gravity waves 



1 INTRODUCTION 

Currently operating or upcoming CMB experiments such 



as EBEX (|Reichborn-Kjennerud et al.||2010|), Spider (|Fil- 



ippini et al.||2Q10|, POLAR BEAR (|The Polarbear Collab 
oration et al.||2010|, Keck (ISheehy et al.||2011| and ABS 



( Essinger -Hileman et al.||2010| will routinely reach the sen- 
sitivity in polarisation required to detect the curl-type pat- 
tern (5-mode) predicted by the simplest models of inflation 
(Dodelso n et al.|2009 ). The predicted amplitude of this sig- 
nal however is comparable or below the predicted signal of 
foreground polarisation over all observationally relevant fre- 
quencies and over most of the sky (Gold et al. 2011). For 
polarisation in particular the foreground signal is dominated 
by synchrotron emission at low frequencies (< 100 GHz) and 
thermal dust emission at high frequencies (> 100 GHz). 

Given the high levels of polarised foregrounds the mis- 
sion planning for these and future experiments requires a 
detailed study of sky coverage to optimise sensitivity to the 
B-mode signal. The impact of the trade off between larger 
sky coverage, depth of observation and the 'cleanliness' of 
or lack of Galactic foregrounds in a patch of sky must all 
be considered. Regardless of how clean the final observed 
sky patch, some level of foreground removal will be required 
for all experiments. Realistic foreground templates based 
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on models and/or observations of the polarisation direction 
and amplitude of foregrounds are very useful when carry- 
ing out this work. However, reliable templates of polarised 
foregrounds at the frequencies relevant to CMB observations 
have been hard to come by and only recently, with Wilkin- 
son Microwave Anisotropy Probe (WMAP) iT-band obser- 
vations have reliable estimates of synchrotron polarisation 
on large angular scales been made. 

Little other polarization data exists at the frequencies of 
interest for future CMB experiments (with frequency bands 
ranging from 90 GHz up to 450 GHz), hence it is necessary 
to model foreground emission by extrapolating the infor- 
mation from existing data. This paper builds on previous 
work presented in FGPolI, which described a model for 
foreground emission due to interstellar dust in the Galaxy. 
The dust model was introduced by |Q'De a (2009) and first 
applied in O'Dea et al. ( 2011 ) for the purpose of studying the 



impact of polarised foreground dust on Spider's ability to 
detect B-mode polarisation. FGPolI gives a detailed expla- 
nation of the dust model and presents a number of full-sky 
template maps at various frequencies. A complete model of 
polarised foreground emission must also include the effect of 
synchrotron emission. This is particularly important for low 
frequency observations i.e. below the CMB 'sweet-spot' at 
100 GHz. This will be the focus of the work reported here. 
Additional components due to spinning dust and free-free 
emission are not thought to give a significant signal in po- 
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larisation and are omitted in our modelling (see for example 
|Macellari et al.| ( |20TT| and |L6pez-Caraballo et al.| ( |20TT| )). 

Synchrotron emission, generated by the gyration of cos- 
mic ray electrons in the Galactic magnetic field (GMF), is 
intrinsically polarised and constitutes the main polarised 
foreground at lower frequencies ( |Page et al. 2007| . While 
the emission from thermal dust is expected to be higher 
than synchrotron emission above around 90 GHz, the sig- 
nal from synchrotron is also not negligible. With the addi- 
tion of synchrotron, this paper provides a complete model 
of polarised microwave foreground emission on large angu- 
lar scales. Here a detailed explanation of the synchrotron 
model is given and full-sky template maps are presented. As 
with previous work the model includes a three dimensional 
description of the Galactic magnetic field on both large and 
small spatial scales. Both polarisation amplitude and angle 
are modelled internally and the templates are scaled such 
that the polarisation amplitude corresponds to a nominal 
value when averaged over the maps. 

This paper is organized as follows. Section [2] briefly re- 
views the mechanism for synchrotron emission. Section [3] 
describes the synchrotron model including the underlying 
GMF, the choice of total intensity template, and the cosmic 
ray density distribution and line-of-sight integration method 
required to evaluate the final Stokes parameters maps. In 
Section [4] examples of template maps produced using this 
model are described and in Section [5] the templates are 
are compared with the WM AP estimated foreground maps. 
Section [6] compares the level of foreground polarisation in 
the templates in coverage areas of a sample of currently 
planned or operating sub-orbital experiments. We conclude 
with some discussion in Section 



2 POLARISED SYNCHROTRON EMISSION 

Diffuse synchrotron emission is one of the dominant Galactic 
foregrounds for CMB observations. The synchrotron radia- 
tion arises when electrons with large relativistic energies are 
accelerated in the GMF. The frequency dependence of syn- 
chrotron emission depends on the energy spectrum of these 
cosmic-ray electrons, as well as the intensity of the GMF 



(see e.g. Rybicki & Light man (1979])). An ensemble of rel- 
ativistic electrons with a power law distribution in energy 
produces a synchrotron emission spectrum that is another 



power law ( Longair (1994)). At GeV energies, where radio 
synchrotron emission peaks, the index of the power law is 
expected to have a range between f3 ~ —3.5 — > —2.5 (Ry- 
bicki fc Lightman|1979 ) in inferred antenna temperature or, 
equivalently, a ~ —1.5 — >> —0.5 in specific intensity. Since 
the spectral index has been seen to vary with position across 
the sky (Hins haw et al.|2007| ), such a power law in antenna 
temperature describing synchrotron emission is only an ap- 
proximation. In addition, the highest energy electrons lose 
energy more quickly resulting in a gradual steepening in the 



power law index at the higher frequencies (Bennett et al. 



2003). 



At microwave frequencies polarised foreground emission 
is dominated by polarised synchrotron and thermal dust 



that are both sensitive to the coherent GMF (Page et al. 
2007). The dominant emission at lower frequencies is po- 



its lower frequency bands provide important constraints on 
polarised synchrotron emission. Synchrotron emission is lin- 
early polarised with direction perpendicular to the projec- 
tion of the GMF on the plane of the sky (see for exam- 



ple Rybicki &; Lightman (1979)). The degree of synchrotron 
polarisation depends greatly on position on the sky and ob- 
serving frequency. Changes in magnetic field direction along 
the line-of-sight leads to a depolarization effect, reducing the 
fractional polarisation degree of synchrotron emission. At 
frequencies lower than ~ 1 GHz, depolarization is significant 
and hence synchrotron polarisation is as low as few tens of 
percent ( jSpoelstra||1984| ) . At CMB frequencies, depolarisa- 
tion is minimal, with the degree of synchrotron polarisation 
being as high as 30 to 50 % in some galactic structure. 

Free-free emission and spinning dust are also thought 
to contribute to the foreground signal in total intensity over 
a range of frequencies. For example, anomalous microwave 
emission around 20GHz has been found in WMAP data, 
with suggestions that this is more likely due to spinning 



dust emission than a flat synchrotron component ( Peel et al. 



2011). However spinning dust and free- free emission are not 



thought to be significantly polarised and their impact on 
final estimates of e.g. the tensor-to-scalar ratio r is minimal 
( Armitage-Caplan et al.|2012| ) at about the 1% level. 

Aside from the spatial dependence of the polarisation 
angle and amplitude, significant uncertainties remain in the 
frequency modelling of synchrotron intensity. For example, 
Ben nett et al.| ( |2003| argue that at higher Galactic lati- 
tudes (in the halo) the spectral index f3 ~ — 3 while in 
the Galactic plane (near star forming regions) f3 ~ —2.5. 
This results in differences in the observed structure between 
WMAP if-band at 23GHz and the Haslam map at 408MHz 



(Haslam et al. 1982) as regions with flatter spectral in- 



dex become more important at higher frequencies. However 
Miville-Deschenes et al.| ( |2008| ) find a lower range of variation 



of the spectral index. This work focuses only on the polar- 
isation fraction and orientation due to the assumed GMF 
model. We assume a simple frequency scaling of the Haslam 
template with a single spectral index multiplying the in- 
ternally modelled polarisation fraction to provide morpho- 
logically realistic templates. A more detailed, possibly pixel 
dependent, frequency rescaling can always be introduced by 
rescaling the template obtained. 

A number of other studies aimed at modelling fore- 
ground emission at microwave frequencies have been carried 
out (see e.g. |Fauvet et all (2010) , |Page et al.] ( |2007| )). An 
extensive study modelling different foregrounds in both in- 
tensity and polarisation over a large range in angular scales 
is the Planck Sky Model (PSM) jDelabrouille et al.|2012| P 



larised synchrotron radiation. WMAP measurements from 



It includes detailed modelling of Galactic diffuse emission, 
including synchrotron and thermal dust emission as well as 
free-free, spinning dust and CO lines. It also includes infor- 
mation on Galactic HII regions, extragalactic radio sources 
and several other sources of emission. 

At present the best template of polarised Galactic syn- 
chrotron emission is that provided by the WMAP if -band 



1 The PSM description was released in the interim betwen FG- 
PolI and this work. The PSM template and model are available 
on a restricted basis and detailed comparisons will be described 
in a future publication. 
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Figure 1. From top to bottom, Q and U Stokes parameter template maps, displayed in Galactic co-ordinates for left: FGPol model of 
synchrotron emission at 23 GHz for the LSA GMF model. These were generated at N s ide = 1024 map resolution based on N^ de = 128 
line-of-sight resolution. The maps are smoothed to 1° and downgraded to N s id e = 64. right: WMAP MCMC synchrotron map for 
comparison. Units are jiK antenna temperature. 



(23GHz) whilst the intensity has been well measured (free 
from CMB contamination) at 408MHz by the Haslam all 
sky survey. As detailed below we use the Haslam maps 
to introduce detailed morphology in our templates since 
the WMAP Maximum Entropy Method (MEM) maps still 
contain a significant noise residual due to CMB contami- 
nation at the smoothing scale adopted (1°). We compare 
the templates obtained here with the WMAP synchrotron 
and dust maps obtained through their Monte Carlo Markov 
Chain (MCM C) best-fit procedur e ( |Gold et al.|20TT| . Simi- 
lar templates ( Fa uvet et al.|2010| have been compared with 
Archeops maps over a limited fraction of the sky at 353 
GHz ( |Benoit et al. 2004) but these maps are not publicly 
available. 



3 MODEL 

3.1 Galactic Magnetic Field Model 

A number of magnetic field models were compared in FG- 
POLl|^| Here we limit the choice to the Logarithmic Spiral 

2 We refer the reader to FGPolI for a detailed discussion of 
GMFs on large and small scales. 



Arm (LSA) model introduced by |Page et al. ( 2007| for use 
in modelling the WMAP data. The model is defined as 



B p = -Bq sin L0o + ipi In 



-) 

PwJ 



cosx, 



B<$> — —Bo cos ( t/jo + ipi In -^— ) cos x , 
B z = Bq sin x , 



(1) 

where p, $ and z are Galacto-centric cylindrical co-ordinates 
with <£, the cylindrical longitude, measured from the direc- 
tion of the Sun, x — Xotanh(z/zo) parametrizes the ampli- 
tude of the z component and zo = 1 kpc. The field amplitude 
is set to .Bo = 3/iG, and we take the distance between the 
Sun and the Galactic center to be 8 kpc. Best-fit parameter 
values obtained by fits to the WMAP K-band field direc- 
tions are ipo = 27 degrees, ipi = 0.9 degrees, and %o = 25 
degrees. The radial scale is set to pw = 8 kpc and the scale 
height is set to zo = 1 kpc. 

Although we focus on large angular scales we also in- 
clude a small-scale random component in our GMF model 
by adding a realisation of a Kolmogorov turbulence field 
with a one-dimensional Kolmogorov energy spectral index 
of —5/3. An injection scale of 100 pc is chosen for the tur- 
bulent realisation with a negligibly small dissipation scale 
compared to the resolution scale. 
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Polarisation Angle WMAP 




Figure 2. Polarisation amplitude (P = \J Q 2 + U 2 ) in antenna temperature (fiK) (top) and angle (bottom) in Galactic co-ordinates for 
the synchrotron emission template (LSA GMF model) at 23 GHz (left) and for WMAP 23GHz (right). The polarisation angle colour 
coding ranges from to 180 degrees. 



3.2 Cosmic Ray Density Distribution 

To model synchrotron emission a three-dimensional model of 
the distribution of cosmic rays in the Galaxy is required. The 
large-scale spatial distribution of the cosmic rays is modelled 
through its density, n cr , and is thought to follow the same 
form as the dust distribution with modified radial and scale 
heights ( |Page et al.|2007| ), 



n cr — n exp 



sech^ 



(2) 



where the height and radial scales are set to z CT = 1 kpc 
and per = 5 kpc. These parameter values were chosen for 
the WMAP analysis of | Page et al.| ( |2007| following work by 
Drimmel fe Spergel| pQQl^ 



3.3 Total Intensity 

Our method aims to predict the polarisation amplitude and 
angle based on a chosen GMF model. The choice of total in- 
tensity for our templates is therefore external to the model 
and can be set to reflect any existing template. For syn- 
chrotron emission we chose to scale the point source cor- 
rected Haslam all-sky survey using a single power law in 
antenna temperature for simplicity. This scaled map is also 
multiplied by the internally modelled polarisation fraction 



template to produce Stokes parameter maps with realistic 
morphology. 

The Haslam template in brightness temperature, with a 
resolution of 0.85°, is scaled to microwave frequencies using 
a spectral index f3 s = —3. Although the map may contain 
residual contamination by free-free emission (see e.g. |Dick-| 
inson, Davies &; Davis ( 2003| ) we assume it is dominated by 
the synchrotron component. The templates can be rescaled 
using any choice of templates in future. 



3.4 Stokes Parameters 

The direction and degree of polarisation from synchrotron 
emission are highly dependent on the Galactic magnetic 
field. To model these we integrate along lines-of-sight us- 
ing the GMF outlined in section |3.1| The full-sky maps 
presented here were obtained using a one-dimensional re- 
alisation of the small-scale turbulent field. When producing 
smaller patches that require much fewer lines-of-sight at a 
given resolution we model the small scale turbulence as a full 
three-dimensional random realisation which preserves the 
spatial correlations implied by the Kolmogorov spectrum. 
The total GMF model is made up of a sum of large-scale 
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(Is) and small-scale (ss) components with 

B r — B r ^ ss -\- B r ^\ s , 
Be = Be, ss + Be, i s , 
B^ = B(j> iSS + B(f)^ s , 



(3) 



where r, and are now Solar-centric spherical polar co- 
ordinates. The polarisation at each point along the line-of- 
sight r is determined by the perpendicular field components, 
Be and B^. 

The Stokes parameters for the synchrotron model are 
then projected out from the three-dimensional model using 
the appropriate line-of-sight integrals, 



/model (0, <j>) = e(u) I 

Jo 

/"^max 

Qmodei(#, </>) = e{v) I n CT (r)p 
Jo 



Urr. 



i(M) 



rr 

<") / 
Jo 



n cr (r)(B (r) 2 + B 9 (r) 2 )dr, 

;g (r) 2 - B 6 (r) 2 )B 2 r 
B 2 



n A r) P 2B ^:^ d r, (4) 



dr , 



where B 2 = B 2 + B 2 + and e is the emissivity as a 
function of frequency, v. As with the dust templates, we 
conform to the default convention applied in the healpd^] 
package (Gorski et al. 2005) regarding the sign of U. 

Having computed the line-of-sight integrals for the 
Stokes parameters we calculate maps of the polarisation di- 
rection, 7, and degree, P, given by 



p(e,4>) 
7 (M) 



V Qmodel + U ? 



model 



/model 

^model 



arctan . 

^ \ ^ model 



(5) 



The final synchrotron template at frequency v is then 
obtained by scaling with the Haslam template 



i sync \ i 
Qsync (^j 
^sync (^j 



^Has(#? 
^Has(#? 
^Has(# 5 



)P((9,( 

|P(0,< 



cos (2 7(0,0)) 
sin (2 7(0,0)) 



(6) 



where /^as is the total intensity of the Haslam map extrap- 
olated to frequency v. 



4 MAPS 

Figure^ shows Q and U Stokes parameter maps at 23GHz 
for the whole sky arising from the model with their am- 
plitudes scaled such that the polarisation amplitude cor- 
responds to that of the WMAP counterpart (also shown) 
when averaged over the maps. The morphology of the po- 
larisation agrees well with the observations with the most 
visible difference being on scales of a few degrees where the 
WMAP estimates are dominated by residual noise. 

The resolution of the healpix maps is Af s id e = 1024 
but the polarisation information is based on a line-of-sight 
integral at an angular resolution of A^ s ^ de = 128, correspond- 
ing to roughly £ ~ 500 in multipole space. We integrate 
along lines-of-sight to the centre of all healpix pixels at a 
given A^ s ^ de from zero out to a maximum distance r ma x of 



: EE : 


: BB : 


Synch (23GHz) 






r>~^Synch (23GHz) \ 




Dust (94GHz)- 







Figure 3. EE (left) and BB (right) fits to the foreground an- 
gular power spectra for the WMAP MCMC synchrotron 23GHz 
and dust 94GHz maps (solid line) compared to the FGPol syn- 
chrotron model and dust model (dashed line). The polarisation 
amplitude of the model templates is matched to the WMAP 
MCMC templates as described in the text. The fit values used 
are given m Table [l] These were calculated from maps masked 
with a union of the WMAP P06 mask and the MCMC flagged 
pixels. 



30,000 pc, with discretisation steps of 0.1 pc. A^ s ^ de is less 
than or equal to iV s id e of the total intensity Haslam map. 

The LSA model is used for the Galactic magnetic field 
model with the same parameters as in FGPolI. The small 
scale field is modeled as Kolmogorov turbulence in ID for 
large patches of sky, with a power spectrum of V(k) oc 
^-(2+3Ar d )/3 w ] iere ]\j d j s the number of spatial dimensions 
of the realisation and k is the magnitude of the wavevec- 
tor. All full-sky templates presented here make use of the 
ID approximation along the line-of-sight for the small scale 
turbulent component of the GMF. Small patch templates 
discussed below are produced with full 3D realisations of 
the field in the (smaller) volume probed by the reduced cov- 
erageQ 

Figure [2] shows maps of P and 7 obtained using this 
choice of resolution and modelling of large and small scales. 
For comparison, maps of P and 7 for the WMAP 23GHz 
MCMC template are also plotted. Differences between the 
templates and observations are mostly due to noise but there 
are also obvious differences in the morphology along the 
galactic plane and around the largest Galactic features such 
as the Galactic centre and North and South Galactic Spurs. 
Some of these differences are related to our choice of to- 
tal intensity template which uses the Haslam maps at 408 
MHz. A comparison between the scaled Haslam map and 
synchrotron templates obtain via the differencing of WMAP 



K and Ka bands were discussed in Gold et aL] ( |2011| ) . Below 
we quantitatively compare the broad features of both syn- 
chrotron and dust full-sky templates with the corresponding 
WMAP MCMC best-fit maps. 



3 See http://healpix.jpl.nasa.gov 



Figure[T]can be compared with Figure 2 in Fauvet et al. (2010) 
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Component 
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WMAP MCMC 














Synchrotron EE 


306 


± 


95 


-0.91 


± 


0.11 


Synchrotron BB 


144 


± 


44 


-0.87 


± 


0.12 


Dust EE 


12.9 


± 


6.4 


-1.06 


± 


0.24 


Dust BB 


6.12 


± 


3.7 


-0.83 


± 


0.29 


FGPol 














Synchrotron EE 


343 


± 


91 


-1.05 


± 


0.09 


Synchrotron BB 


110 


± 


29 


-0.73 


± 


0.08 


Dust EE 


1.38 


± 


0.33 


-0.12 


± 


0.07 


Dust BB 


1.70 


± 


0.37 


-0.22 


± 


0.06 



Table 1. Foreground power law + white noise fits of WMAP 
MCMC and FGPol template spectra outside the combination of 
P06 mask and MCMC flagged pixels. There is good agreement in 
both EE and BB for synchrotron between the FGPol template 
and the WMAP MCMC synchrotron component map. The FG- 
Pol dust template shows a significantly shallower spectrum than 
the WMAP MCMC component map indicating relatively more 
structure at large angular scales. 



5 COMPARISON WITH WMAP TEMPLATES 

The WMAP satellite observations provide full-sky maps of 
temperature and polarisati on in five frequenc y bands be- 
tween 23GHz and 94GHzp]( |jarosik et al.|201l| ). The polari- 
sation maps contain important information on Galactic fore- 
ground emission and hence provide an important test of our 
model of Galactic synchrotron radiation. Synchrotron radi- 
ation dominates the measured signal in the lower frequency 
bands, and we use the best fit WMAP 23GHz synchrotron 
templates generated by an MCMC fit for comparison with 
our model. Thermal dust emission dominates higher fre- 
quency bands, and for comparison with our dust model we 
use the 94GHz dust templates generated from the MCMC 
best fit values. 

We use the WMAP MCMC 'base' fit which includes 
three power law foregrounds: dust, synchrotron and free-free 
emission as well as a contribution from CMB. These maps 
are smoothed at a scale of 1° and have been downgraded to 
N side = 64 before the MCMC fit is carried out. The WMAP 
team performed a combined MCMC fit to their five bands 
at a resolution of N s id e = 64. The pixel noise is calculated 
at A^ s ide = 512 and downgraded to N s id e = 64. 

FGPol templates are normalised to the polarisation 
amplitude P = y/Q 2 + U 2 of the full sky WMAP MCMC 
template with no masking or smoothing. The model maps 
are generated at N s - 1( \e = 1024 for the large scale resolution 
and N^ ide = 128 for the small scale line-of-sight resolution, 
then smoothed with a Gaussian beam of 1° and degraded to 

iVside = 64. 

Angular power spectra Cf x for XX = TT, EE, and 
BB of both templates were calculated after masking with 
the P06 polarisation mask (Pag e et al.|2007 ) combined with 
a mask of pixels nagged by the WMAP MCMC process. 
The spectra are corrected for sky fraction / s k y effects and 
for their respective pixel and beam smoothing functions and 
then fit for a power law in multipole £ with an additional 




Dust FGPOL xWMAP (94GHz) 

— Q 

--- u 



20 40 

Galactic cut 0„ 



Figure 4. Comparison of the correlation coefficient 7Z between 
the FGPol maps and WMAP best-fit templates for both syn- 
chrotron and dust. In order to look at the correlation on large 
scales we smooth these maps to a resolution of 10°. The cor- 
relation is calculated by including areas outside a range of cuts 
in Galactic latitude. The trend shows that both dust and syn- 
chrotron models are highly correlated with the best-fit maps at 
high Galactic latitudes whilst the correlation falls rapidly at low 
galactic latitudes for the dust comparison. The WMAP best- 
fit dust map however has an even larger noise residual than the 
synchrotron map and this is expected to reduce the correlation 
significantly. 



white noise component 
£U + 1 



2tt 



d = A£ m + £(£ + 1)AT 



where A is the amplitude of the foreground component, m 
is the index and N is the noise amplitude. This procedure is 



similar to that carried out by Gold et al. (2011). Although 



http:/ /lambda, gsfc.nasa.gov 



the scatter at large angular scales is non-Gaussian and some- 
what correlated by the sky cut we adopt a very simple as- 
sumption for sample variance in the power spectra by dis- 
regarding correlations between multipoles and assuming a 
Gaussian scatter given by the sample variance for each C xx . 
The WMAP derived fits also make use of diagonal Fisher 
error values. 

The fits allow us to quantify the scaling of the angular 
power spectrum for both templates as a function of multi- 
pole £ whilst allowing for any residual noise and/or pixeli- 
sation effects. They can also be used as a quick guide for 
the level of foreground contamination at different frequen- 
cies on large angular scales either on the full sky or on small 
patches. Figure [3] shows the resulting power law fits in both 
Cf E and Cf B for the P06 masked FGPol and WMAP 
MCMC synchrotron templates at a frequency of 23 GHz. 
Also included are the results of the same procedure applied 
to the dust FGPol and WMAP MCMC templates at 94 
GHz. 

The fit values, excluding the noise amplitudes, can be 
found in Table ^ The synchrotron templates agree well with 
the WMAP MCMC maps in both amplitude and angular 
dependence whereas there are significant differences between 
the FGPol dust template and the WMAP MCMC map at 
94 GHz. 

We also attempt to quantify the level of correlation be- 
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Figure 5. Comparison of the scatter between the FGPol maps 
and WMAP best-fit templates for synchrotron. This is done for 
the full-sky templates including pixels outside the combination 
of the P06 polarisation mask and WMAP MCMC flagged pixels. 
The dust correlation density is omitted as it is dominated by the 
variance in the WMAP map due to residual noise. 



tween the WMAP MCMC maps and FGPol templates. We 
do this using two separate measures. The first analyses the 
level of pixel-to-pixel correlation between maps calculated 
for the area of the sky outside a given galactic latitude cut. 
The correlation coefficient 1Z is given by 



n(o g ) = 



E P (w P -w P )(^ P 



(8) 



where W and T are the /, Q, or U Stokes values of the 
WMAP and FGPol maps respectively and the index p sums 
over all pixels outside the cut at latitude ±0 g . The result, 
for both dust and synchrotron Q and U Stokes parameters 
is shown in Figure [4] The analysis shows that both dust 
and synchrotron templates are highly correlated with the 
WMAP best-fit foreground templates at high Galactic lati- 
tudes. Whilst this is also true for synchrotron at low Galac- 
tic latitudes, the dust model fails to reproduce the observed 
morphology well at latitudes below ~ 30 degrees. This is 
not surprising since thermal emission by dust particles is 
more susceptible to the detailed structure in the Galactic 
disk with even large angular scales being influenced by tur- 
bulence and/or existence of individual clouds. 

We also look at the correlation in terms of scatter of the 
pixel values in Q and U Stokes parameters for the WMAP 
MCMC maps versus the FGPol templates in both dust and 
synchrotron. All pixels outside the P06 mask and MCMC 
flagged pixels are included and the scatter density is shown 
in Figure [5] as two contours encompassing 68% and 95% 
of pixels. We only show the synchrotron correlation density 
since the dust one is found to be dominated by the larger 
WMAP variance due to residual noise and is not informa- 
tive. 



6 FOREGROUND AMPLITUDES IN 
SUB-ORBITAL SKY PATCHES 

We also examine the amplitude of foreground contamination 
in smaller sky areas being targeted by a sample of three cur- 
rently operating or planned sub-orbital experiments; EBEX 



, Spider ( 
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rays which observe the same field. Angular resolution and 
sensitivity for the three experiments are varied but they are 
all targeting the detection of BB power either on large an- 
gular scales that are free of lensing effects or, as in the case 
of EBEX, on smaller angular scales where the lensing effect 
dominates the BB signal. 

We generate high resolution templates with full three- 
dimensional modelling of the turbulent small-scale GMF 
over the regions of expected coverage for the three experi- 
ments. For the dust templates, the Q and U components are 
normalised so that the average polarisation fraction outside 
the area defined by the WMAP P06 mask is 3.6%. The cov- 
erage areas are outlined in Figure [5] with Spider targeting 
the largest area with / s k y ~ 0.1, EBEX targeting the small- 
est patch contained in the Spider area with / s k y ~ 0.01 
and Bicep2/Keck targeting the southern most patch with 
/ sky - 0.03. 

In order to compare the relative contamination by fore- 
grounds in relation to the relevant signal we analyse the an- 
gular power spectrum Cf B for both dust and synchrotron. 
Due to the uncertainties involved in predicting the small 
scale signal we focus on large scales only with templates 
smoothed to a common resolution of 1° and rely on ex- 
trapolating a power law to scales larger than £ ~ 200 in 
comparing with the expected signal. 

The analysis on small areas of the sky such as these is 
complicated by the significant correlation induced by the cut 
on spherical harmonic coefficients. The high level of correla- 
tion would result in significant biases if the same power law 
fitting procedure as used in Section [5] were carried out. To 
avoid this problem we estimate the overall amplitude of fore- 
ground contamination by averaging in pixel space assuming 
a fixed power law in £ corresponding to our previous near 
full-sky analysis. 

In practice we calculate the variance in both Q and U 
for each patch and assume a relation between the variance 
and angular power spectrum of the form 



47T 



£(2^ + i)c,b; 



(9) 



with the signal angular power spectrum modelled as Ct = 
A£ m in accordance with |7| and with index m set to the 
corresponding near full-sky best-fit value (see Table [lj. We 
take tax = 128 and model the smoothing Bg applied to the 
templates as a Gaussian beam with FWHM 1° multiplied 
by the pixel window function at the working HealPix reso- 
lution A^ s ide = 64. We then invert the relation ([9]) to obtain 
an 'average' polarisation angular power spectrum amplitude 
A, effectively assuming that power is equally distributed be- 
tween EE and BB. 

The results are summarised in Figure [7] for a single ref- 
erence frequency of 150 GHz as this is being included as 
an observing frequency in all experiments being considered. 
The model power spectra for each patch are shown in ther- 
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Figure 6. Comparison of patches targeted by a variety of sub- 
orbital experiments. Spider targets the largest area (/ s k y ~ 0.1), 
EBEX targets the smallest patch (/ s k y ~ 0.01) and Bicep2/Keck 
targets the southern most patch (/ s k y ~ 0.03). 



modynamic temperature in order to compare directly with 
the expected BB signal for a tensor-to-scalar ratio r — 0.1. 
Both primordial and lensing contribution to the BB signal 
are shown. 

The amplitude of foreground contamination varies by 
roughly an order of magnitude between the area targeted 
by different experiments. In particular the area targeted by 
EBEX seems to be very clean with the foreground signal 
reduced by an order of magnitude compared to the areas 
targeted by Spider and Bicep2/Keck. This agrees visually 
with the impression given in Figure [6] 



7 CONCLUSIONS 

We have presented templates for polarised emission from 
synchrotron radiation within our Galaxy using a 3D model 
of the Galactic magnetic field and cosmic ray density dis- 
tribution. From this model, maps of polarisation amplitude 
and angle are calculated which are then combined with to- 
tal intensity measurements from the Haslam 408MHz all-sky 
radio continuum survey to provide template maps. 

We have compared the FGPol templates obtained from 
this model with data from the WMAP satellite for both 
synchrotron and dust emission. We find that the synchrotron 
template agrees qualitatively with the observations whereas 
comparison of the dust template is complicated by the large 
residuals present in the WMAP estimates. 

We have also looked at foreground contamination levels 
in patches that will be targeted by upcoming experiments 
and found that our model predicts significant differences of 
up to an order of magnitude in the foreground contami- 
nation of different patches. The level of contamination will 
dominate the ability of various experiments to achieve their 
target sensitivity with respect to the B-mode signal being 
searched for. 

As more polarisation data becomes available the com- 
parison between the model and observations will become 
more quantitatively precise. In particular future Planck data 
releases will provide high signal-to-noise Q and U maps at a 
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Figure 7. Foreground amplitude for synchrotron (left panel) and 
dust (right panel) calculated from patches targeted by various 
suborbital experiments compared with theoretical EE and BB 
spectra for r = 0.1. The index of the power law used is from 
the corresponding near full sky power spectrum fits. Also shown 
are the near full sky best fit spectra for the FGPOL synchrotron 
template from Table ^ along with best fit spectra to the dust 
templates with Q and U normalised so that the polarisation frac- 
tion is 3.6% outside the WMAP P06 mask. The amplitudes were 
calculated from our FGPOL dust and synchrotron templates at 
150GHz calculated from maps generated at iV s ide = 1024 for the 
large scale resolution and iV^de = 128 for the small scale line-of- 
sight resolution, which are then smoothed to 1° and downgraded 
to Af s ide = 64. Units are /ii^cMB- 



number of frequencies and we will be able to refine our model 
based on them. Indeed, in future, it should be possible to 
learn much about the Galactic magnetic field itself by fitting 
the (many) model parameters to actual data. This will shed 
light on many aspects of our Galaxy's physical model that 
are still poorly understood. 
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